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Bose-Einstein Condensation in Competitive Processes 
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We introduce an irreversible discrete multiplicative process that undergoes Bose-Einstein conden- 
sation as a generic model of competition. New players with different abilities successively join the 
game and compete for limited resources. A player's future gain is proportional to its ability and 
its current gain. The theory provides three principles for this type of competition: competitive 
exclusion, punctuated equilibria, and a critical condition for the distribution of the players' abili- 
ties necessary for the dominance and the evolution. We apply this theory to genetics, ecology and 
economy. 



PACS numbers: 03.75.Nt, 87.23.-n, 87.23.Kg, 



1.65, 



Competition occurs when two or more players such as 
organisms, individuals or companies strive for common 
but limited resources. It plays a significant role in biolog- 
ical and social activities, and is the basic idea underlying 
evolution. Recently Bianconi & Barabasi demonstrated 
that their competitive complex network model behaves 
analogously to a Bose gas 1]. A few nodes recruited 
many connections and acted like hubs, while most others 
had sparse links. This phenomenon was described in the 
mathematical framework of Bose-Einstein condensation 
(BEC). Here, we analyze a discrete multiplicative process 
that also exhibits BEC. This generic competition model 
is defined by the following conditions at each time step: 
1. A player's gain is proportional to its ability and to its 
gain at the previous time step. 2. A new player joins the 
competition. 3. Players compete for a fixed total amount 
of resources. The constant addition of players is the basis 
of the model's open and irreversible characteristics and 
makes the process suitable for the analysis of natural and 
social competition. We show that the cumulative gain of 
each player is related to the Bose distribution, except 
for that of two players. One is the pioneer. The other 
is a high-ability player that dominates most resources 
under a certain statistical condition for the distribution 
of the players' ability. We investigated the implication of 
BEC in this mathematical game for different competition 
models. Applied to the competition among clonal strains 
of evolving asexual bacteria, BEC in an irreversible pro- 
cess provides theoretical grounds for the principle of com- 
petitive exclusion and the hypothesis of punctuated 
equilibria 0- Furthermore, the existence of a transition 
point elucidates the condition necessary for evolutionary 
development, and explains the observation that a high 
mutation rate is not necessarily advantageous for adap- 
tive evolution 0, 0|. Finally, we investigate possible ap- 
plications to ecological and economic activities. 



Before we present our competition model, we introduce 
a novel, general framework for how a function ip(e) for- 
mally maps into a Bose system. Define a function ip(e) 
that satisfies the following conditions for an arbitrary C 2 



density function g (e) (> 0) defined on e e [0, e max ], 



deg(e) 



deg{e)ip(e) = m + yo, 

deg{e)e- 0e tp{e) = M. 



(1) 
(2) 
(3) 



The r.h.s. of these equations and j3 in J3J| are all positive 
constants. Then ip (e) is given by 



<p(e) = 



(4) 



where e Q = m/M. For, dividing J5| by yo, then sub- 
tracting (JIJ yields 



deg(e) = N, 



yo 



(5) 



where N = m/y - From ((3J xN/M— (J5J and the funda- 
mental lemma of the calculus of variation, we obtain 
Hence, the so-called occupation number in JSJ, 



n(e) 



<P (e) ~ Vo 

yo 



o0e+a _ I ' 



(6) 



becomes the Bose distribution, or, equivalently, the sys- 
tem tp (e) maps into a Bose system. From (@J, such a 
system may be obtained from the sum of a geometric 
progression with ratio e~^ £ ~ Q < 1. This motivates anal- 
ysis of the following competition process. 
Consider a simple multiplicative process 



z/i(t + 1) = aib(t)yi(t), 



(7) 



where yi(t) is the gain of the ith player at time t and 
di is its ability, a positive and time-independent random 
variable chosen from a distribution p(a). From J2J, a 
player's gain is proportional to its ability and the gain of 
the previous time step. Furthermore, a new player joins 
the race at every time step with initial value yi(i) = yo, 
which means we consider t + 1 simultaneous equations at 



2 



time t. The only exception is the first player (pioneer), 
whose initial value is yo+m. Finally, the population races 
for a fixed amount of resources, modelled by the normal- 
ization b(t) in J7J which is defined by b(t) = m/M t , where 
M t is given by 



M, 



(8) 



This assures that the total gain distributed among t play- 
ers at time t is m for t > 0. As a new player joins the 
race with yo, the total gain of the population is fixed to 
m + Ho every time step. We now consider the time evo- 
lution of players except for the pioneer. The gain of the 
ith player at time t is given by yi(i) — f2j(i)yoj where 



't-i 



fh(t) 



l[b(t') 



(9) 



Without loss of generality, we can set di < 1 (discussed 
later), and we let a, = e~@ ei , thus introducing a for- 
mal parametrization of the ability distribution by the 
inverse temperature /3(= 1/T). The random variable 
e € [0,e max ] is chosen from g(e), a state density func- 
tion which defines the system. From (JjJ, we have 



In fii = (t-i){-0ei - (ln6». 



(10) 



using lnn*r = >(i') - D^- In 6(f) ~ (i-i) (In 6). By 
assuming stationarity of b(t), discussed below, we let its 
time average be 



(In 6) 



(11) 



The normalization factor a keeps the total gain constant 
throughout the dynamical process. Its role is analogous 
to that of the chemical potential of a quantum gas in- 
troduced for the conservation of particle number. From 
JHJ and dJ, we obtain yi {t) = e^^'^-Vyo. Given 
e -i3 ei -a < t]le cumu i a ti V e gain <p t (e<) = Yl\'=iVi (*')> 
converges to tp (ej) in Q. Hence, the normalization of 
¥>* (fi) by the initial value yo, 



n t (ej) 



<Pt (ex) - 2/o 



(12) 



takes the form of the Bose distribution in the thermo- 
dynamic limit (t — > oo). Indeed, in this limit, M t obeys 



/oo 
deg(e) Vc^(/). (13) 
* — * 

t'=0 



Since the stationarity of b (t) yields lim In b (t) — (In b) , 

t — >oo 

we have lim m/M t = e~ a from 111|) . By substituting 

£ — >oc 

this into l|13|l . we obtain an equality, 



f deg (e) -5— ^ = N, 



(14) 



where iV = m/j/o- The dynamical system JJJ thus for- 
mally maps into a Bose system in the thermodynamic 
limit. 

Comparing this system to a quantum gas, adding the 
gain m at each step corresponds to adding m Bose par- 
ticles. This gain is then distributed among players with 
different abilities; these abilities correspond to energy lev- 
els. Different from a quantum gas, a new energy level is 
introduced at every step (with initial occupation of yo), 
the gains (particles) do not move to different energy levels 
once they are allocated, and they are not quantized. De- 
spite these differences, the system preserves a basic fea- 
ture of a Bose gas: gains tend to accumulate where other 
gains are already present. Thus, condensation analogous 
to BEC is expected at low temperatures. 

To study this prediction, we simulated the multiplica- 
tive process (JTJl, adopting a standard density function, 



9(e) = C a 



(15) 



where C„ = er/e^ ax (a > 1). We found two distinct 
phases for the distribution of normalized cumulative gain 
(JEJ- At high temperature, n t (ej) obeys the Bose dis- 
tribution (FlC^-a). Our assumption of stationarity of 
b {t) is confirmed numerically since nt (ej) agrees well with 
the Bose distribution with a being provided on the as- 
sumption of stationarity. At low temperatures, the co- 
ordinated distribution breaks down: one player with low 
energy (but not necessarily the lowest energy because of 
the normalization) dominates a large fraction of the cu- 
mulative resources. Note that the stationarity of b(t) no 
longer holds in this domain. This Bose-Einstein conden- 
sate exists only below the critical point a = 0. For, if 
a is negative, f[ei) of the player with ei < —a/ (3 does 
not converge. Its normalized cumulative gain cannot be 
incorporated in the integral in (|14Jl : rather it has to be 
added as an extra term. The critical temperature for this 
transition can be predicted in the usual way. Change of 
variables y = (3e and use of a = in 1(14^1 yields 



N 



(/^c^max) Jo 



eJJ - 1 



dy. 



(16) 



By approximating the upper limit of integral in l|l(il) by 
infinity, the critical temperature T c (= 1//3 C ) is given by 



T (i) 



{N-'T (a + (a)} - 



(17) 



for N <C 1. When N 3> 1, this approximation is not valid 
because of the high critical temperature. For N ^> 1, one 
uses e v ~ 1 + y in (|16[) and obtains 

Ti 2) ~ e max (1 - cj- 1 ) N. (18) 

The transition can be easily seen in the occupation by 
the most capable player, defined by tp t (e ro in) j J2j ft ( e i) > 
where e m ; n = min {e^}. As the temperature decreases be- 
low T c , its occupation dramatically increases, supporting 
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Temperature (T - 1/p) 



FIG. 1: (a) Normalized cumulative gain (II 12 P for m — 5 and 
T < T c (left), T « T c (center) and T > T c (right). Solid line 
is the Bose distribution with a from H1H . The plus sign indi- 
cates a pioneer. Plots of the last ten entrants were excluded, 
(b) Numerical calculation of |a| with a from 1111 . averaged 
over 500 trials. Symbols +, o, x indicate m — 5, 10 and 20. 
Change in the sign of a indicates the transition. Triangles 
on the abscissa are the transition temperatures from 1181 : 
T c « 3.33, 6.67 and 13.3. Shown in inset is the dependence of 
Tc on exponent a for m = 20 (dots) and the analytical solution 
from l|18|l . solid line, (c) Cumulative occupation by the most 

capable player was calculated by ipt (e m in) / Ylj ft ( e j) > where 
e m i n = minjej} and t = 200. Symbols as in (c). Dashed line, 
occupation by pioneer for m — 20. In all simulations, yo = 1, 
cr = 3 and e max = 1. 



the prediction (FIG. QJ-c). The time evolution of each 
player's gain also shows characteristic differences above 
and below the predicted T c , At T > T c , the gain of each 
player monotonically decreases. Below T c , time evolution 
is not monotonic (FIG.|3-b and area in FIG. (2J-c). This 
is particulary interesting because this BEC occurs in an 
irreversible process. In contrast to BEC in a quantum 
gas in equilibrium, the condensation state dose not im- 
ply the dominance of resources by a single player at the 
lowest energy throughout the whole competition process. 
Instead, there are replacements of dominant player upon 
the entrance of a player with lower energy level. 

We verified the existence of a BEC analogue in a 
discrete multiplicative process, as was shown in a con- 
tinuous model 1]. To see its applicability to real- 
world competition, we mention that Q is a natural ex- 
tension of Malthusian population increase: population, 
when unchecked, increases in a geometrical ratio, that is 
Hi(t + l) =a,iy (i), where 

®i ^ [^min ; ^max 

] and a min > 1. 



We introduced normalization to model limited resources, 
obtaining Q with aj replaced by Sj. Note that we can 
substitute 5,; = a max ai (i.e. < 1) without loss of gen- 
erality because a max is common in the nominator and 
the denominator in J7J). This enables us to apply @ to 
population dynamics and even economic activities. 

Application to Genetics: We consider competition of 
clonal strains of asexual Escherichia coli serially prop- 
agated on glucose-limited medium. For this aim, we 
slightly generalize model (7J such that we allow I (> 1) 
new players to join the competition with initial value yo 
at each step. From a similar argument, we obtain 1141 . 
where N = m /yol . Assigning yo — 1, we have I mutants 
(with usually different abilities) at each time step • We 
chose this procedure because it is unlikely to have yo > 1 
mutants of the same ability. The original m + yol wild 
types propagate m wild types and produce yol mutants of 
novel abilities. Then the mixed population of m + yol is 
normalized to m, and generates yol mutants (FIG.[21-a, c). 
Eventually all wild types are swept away and replaced by 
mutants. This is much faster in BEC than at high tem- 
peratures (see reduction of pioneer occupation in BEC in 
FIG. nj-c). Also at low temperatures inferior strains are 
quickly expelled by a superior strain, which is analogous 
to the competitive exclusion principle 0, @ • Subsequent 
transient replacement of the dominant mutants supports 
the hypothesis of punctuated equilibria Q- Weighted 
mean ability (M t in JSJ if I = 1) shows step-like evolution 
at low temperatures (FIG.|3-d), but this effect diminishes 
with increasing temperature. Now, new mutants all have 
similar abilities, competition is crowded and severe, and 
the sudden upsurges of the weighted mean ability of the 




FIG. 2: (a) Schematic diagram of propagation of E. coli. 
Symbols, o: wild type (WT), others: mutants, (b) State den- 
sity l|15|l (dashed line) with a = 3 and e max = 1 used in 
simulation and corresponding ability distribution (solid line) 
p (a) = g (e) l^l for T = 1 /3, 1 and 3 from left to right. Note 



a 6 e 



-Per, 



.1]. 



Time evolution at T = 1 (m = 200, / = 2, 
yo — 1). As T c ~ 66.6, the process is in BEC. Graded areas 
indicates population of each strains. WT abilities were set 
e = 0.5, which is slightly lower than (e) w 0.76 for a — 3 since 
WT is likely to fit to environment better than chance, (d) 
Weighted mean ability shows punctuated equilibria at T = 0.1 
(bottom) and T = 1 (center). At T = 100 > T c (top), sta- 
tionarity (no evolution) is observed. 
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FIG. 3: (a) Market share of all nine different razors from The 
Gillette Company introduced over the last 3.5 decades Q- 
(b) Simulated market share from J7J, using yo = l,m = 9. 
Rather than assuming a specific state density, we assigned 
a = 0.025,0.04, 0.043 ... etc. for successive products, with 
each newly introduced product having higher ability than the 
previous one. The specific values were selected to best fit the 
empirical data. In years without a product introduction by 
Gillette, the simulation introduced "products" with randomly 
drawn lower abilities (a <C 0.025) that are not competitive 
and disappear (traces adjacent to abscissa). Share of the jth 
product at time t was defined by yj{t) /(m + yo) x 45 as 45% 
is the mean total share in (a). We considered the Carbon 
razor as the pioneer because it was the dominant product at 
the beginning of the period for which we have data. The 
introduction year of "Prem Double Edge" was not available, 
thus we assigned 1963 arbitrarily to fit the plot. 



population are smeared out. Finally, at sufficiently high 
temperatures, no mutant becomes an absolute winner; 
evolution ceases when the population lacks prominent 
players. Note that the stationarity of weighted mean 
ability confirms the stationarity of b(t), which in turn 
allows the argument below Qll[). 

A striking statement of the theory is the prediction of a 
singular point on the emergence of evolution. Most muta- 
tions are likely to be deleterious , which was simulated 
by the use of an increasing state density i|15fl so that infe- 
rior mutants are likely to occur (FIG.[2-c). In this case, 
the critical temperature is predicted by i|18f) . Through 
its dependence on N, T c is found to be dependent on the 
mutation rate yol/(m + yol) = (N + I) 1 . Therefore, at 
a given temperature, decreasing mutation rate generates 
evolutionary development. A sufficiently low mutation 



rate not only prevents deleterious offspring but is neces- 
sary for adaptive evolution itself. 

Although introducing the notion of temperature and 
the state density of a quantum gas was instructive, our 
model does not require these constructs. Our conclusions 
hold as long as a sharply decreasing ability distribution 
p\{a) exists on a 6 [a m j n ,l] where a m - m > 0, with A a 
parameter, such that the integral l|T3J). 



/(A, a) 



dapx (a) 



1 



N, 



(19) 



L e Q - 1 

is finite. Since / (A, a) < I (A, 0) = N c (A), we have a fi- 
nite N c and we can determine the corresponding critical 
mutation rate for a given A. Since iV c is finite, the process 
is in condensation for N > N c , while evolution ceases for 
N < N c . In those cases in which the critical mutation 
rate is in the physiological range, this argument explains 
"clonal interference" [8j which is an empirically observed 
constraint on the role of mutator alleles that rises to high, 
frequency in association with beneficial mutations 0, 0| • 
The critical point, however, only determines the occur- 
rence of adaptive evolution at t — ► oo, and a low mu- 
tation rate slows the rising of dominant mutants. This 
implies the existence of an optimal mutation rate for a 
transiently changing environment. 

Application to Ecology and Economy: In ecology, the 
Lotka-Volterra equation predicts that coexistence be- 
tween strongly competing species is impossible, a result 
confirmed in laboratory experiments on Paramecium Q 
though applicability in complex ecosystems is uncertain 
. These results established the principle of competitive 
exclusion] or "one species, one niche." Our process mod- 
els interspecific competition with completely overlapping 
niches, and BEC is a natural extension of the competi- 
tive exclusion principle in an irreversible process. This 
principle may also apply to economic activities if a prod- 
uct or a company satisfies the ecological conditions of the 
exclusion principle [Toj . For instance, products sold for 
a specialized market may compete within their niche in 
the marketplace. Figure 0-a shows the market share of 
different razors successively introduced by The Gillette 
Company, illustrating a much-studied case of new prod- 
ucts eroding the market share of their own predecessors 
j6f . Aspects of the dynamics of product replacement may 
thus be captured by our competitive process (Fig. |3J-b). 
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